Prom quantum circuits to adiabatic algorithms 
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This paper explores several aspects of the adiabatic quantum computation model. We first show 
a way that directly maps any arbitrary circuit in the standard quantum computing model to an adi- 
abatic algorithm of the same depth. Specifically, we look for a smooth time-dependent Hamiltonian 
whose unique ground state slowly changes from the initial state of the circuit to its final state. Since 
this construction requires in general an n-local Hamiltonian, we will study whether approximation 
is possible using previous results on ground state entanglement and perturbation theory. Finally we 
will point out how the adiabatic model can be relaxed in various ways to allow for 2-local partially 
adiabatic algorithms as well as 2-local holonomic quantum algorithms. 

PACS numbers: 03.67.Lx 



1. INTRODUCTION 



Adiabatic evolution as a quantum computation model 
has attracted much attention since its introduction by 
Farhi et al The basic idea is the following: Start 
with a Hamiltonian whose ground state is easily reach- 
able and prepare our state in the ground state. Change 
it slowly to a new Hamiltonian that encodes the solution 
of the problem and maintain a large energy gap between 
the ground state and the excited state that the evolving 
state couples to. The Adiabatic Theorem then guar- 
antees that the resulted state will be very close to the 
ground state of the new Hamiltonian. The original form 
of the Hamiltonian considered in [l| is a straight-line in- 
terpolation: H(s) = (1 - s) H initial +sH 'final- Recently, it 
was proved that any standard quantum circuit, specified 
by a sequence of unitary operators, can be implemented 
as an adiabatic evolution of this form0,0|. The authors 
use computational complexity techniques developed for 
proving the QMA-completeness of the k-local Hamilto- 
nian problem (Kempe et al Q achieved the case for k=2); 
the evolving state encodes the entire computational his- 
tory. Roughly speaking, they construct a Hamiltonian 
whose ground state is the superposition of all the stages 
in a given circuit. If the circuit has depth L, the time 
required to obtain this ground state is 0(1/ L 6 ) for the 
3-local Hamiltonian and there is a 0(1/ L) probability of 
obtaining the final state of the circuit given this super- 
position. On a seemingly unrelated note, Farhi et al 
showed after Q that if we do not restrict adiabatic evolu- 
tion to the straight-line path and add terms that vanish 
at the endpoints, we may be able to turn an inefficient 
computation into an efficient one. A general method for 
finding an efficient path is however not known. In light 
of these two developments, we may ask - Can we always 
find an efficient adiabatic evolution path, not necessarily 
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of the straight-line fom, for problems efficiently solvable 
by quantum circuits such that we directly obtain the de- 
sired final state? Starting with this question, we will 
present several variations on how to implement the stan- 
dard circuit model by adiabatic evolution. 

The first result we show is that once we specify a) the 
unitary transformation that takes the eigenstates at the 
beginning to those at the end of the evolution, and b) 
how we want the eigenvalues to evolve, we can imme- 
diately derive the form of time-dependent Hamiltonian 
required without the use of any ancilla qubits. Con- 
ceptually the simplest example is a time-dependent sim- 
ilarity transform, U(t)Hi n itiaiU(t)^ . This observation 
should allow us to engineer Hamiltonians according to 
computational needs. However, a Hamiltonian of the 
type ...U2U\HinitiaiU\U\... can be highly non-local even 
if £/i, t/a... and ilinitiai have simple local forms. Since we 
are interested in the ground state, we should ask whether 
it is possible to find an local approximation. 

It turns out that while approximations are possible to a 
certain extent, there is much constraint. We will demon- 
strate this point in two steps. First we will make use 
of the results by Haselgrove et alp, which show how 
the entanglement of the eigenstates of a Hamiltonian is 
related to what bodies in the system each term in the 
Hamiltonian acts non-trivially on. Intuitively speaking, 
if an eigenstate shows strong correlation between bodies 
which the Hamiltonian does not directly couple, i.e. act 
nontrivially on all as a tensor product, the Hamiltonian 
cannot distinguish very well between such a state and 
other similarly entangled states that are orthogonal to 
it. This results in a small energy gap. Since a quantum 
circuit can generate highly correlated states, when we 
want to make them ground states of a Hamiltonian in an 
adiabatic algorithm, they will be difficult to approximate. 
Then, as an explicit example, we will use the approxima- 
tion method developed in derived from perturbation 
theory, and apply it on our construction. We will see 
that we could indeed make a local approximation under 
the constraint implied by 6], but the resulting evolution 
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can be inefficient. 

Next we iook at how this approach of transforming the 
Hamiltonian adiabatically is related to the manipulation 
of geometric phase. We start by asking, given local ap- 
proximation is difficult: Why is the adiabatic model more 
demanding than the basic circuit model, for which 2-local 
Hamiltonians easily suffice with U = exp(iH2-i OC ait) for 
each gate? There are at least two crucial differences be- 
tween the two models. 

I) The adiabatic model keeps track of exactly where 
the state is at every moment throughout the 
evolution and penalizes any deviation, while the 
circuit model keeps no information about the state 
at all. This makes the former more resistant to 
error. 

II) The adiabatic model allows time variability in the 
application of the Hamiltonian, while the basic cir- 
cuit model requires precise pulse timing. 

The word "adiabatic" itself only suggests property II) 
above, so if we are willing to relax property I), we would 
have much more freedom to design our Hamiltonian. The 
main issue we need to deal with, as we will show, is 
the geometric phase. Suppose we implement the Hamil- 
tonian U(s)Hi n iti a iU(sy without making sure that we 
start with the ground state or any eigenstate of the ini- 
tial Hamiltonian. Instead of having U (s) applied on 
our initial state, there will be further transformation 
due to the relative phases accumulated between different 
eigenstates. The dynamical component of the phase is 
straightforward to cancel out, but the geometric compo- 
nent is more subtle to calculate. To avoid having to can- 
cel the geometric component, we may either make sure 
we always start with an eigenstate of the Hamiltonian, 
or take advantage of the geometric phase to implement 
the desired transformation. This would naturally lead us 
to the holonomic quantum computing model developed 
by Zanardi et alQ, which in fact precedes the adiabatic 
algorithm of []|. 

Let us give a lightning review of the idea of holonomic 
quantum computing (HQC). Wilzcek and Zee introduced 
in H the observation that if a Hamiltonian with degen- 
erate eigenstates goes through a cycle adiabatically with- 
out changing the degeneracy of each level, the degenerate 
subspace can be viewed as a gauge group on the mani- 
fold corresponding to the parameter space of the Hamil- 
tonian. After each cyclic evolution, an arbitrary state 
in the degenerate space will undergo an unitary trans- 
formation depending on the path taken; the set of all 
possible such unitary transformation given a parameter 
space that specifies the Hamiltonian is called the holon- 
omy group, and the parameter space is often called the 
control manifold. Elements of the group generally do not 
commute, so the transformation is called the non-Abelian 



geometric phase. Zanardi et al ,7| applied this idea on 
quantum computing by choosing initial Hamiltonians for 
which the computational states are completely degener- 
ate. Transformations are then applied by holonomy. In 
addition to time variability the geometric nature (such 
as dependency on the area of the loop) also gives HQC 
some resistance to errors. 

Seeing HQC as a generalized adiabatic model brings us 
many new insights. First, we can apply the local approx- 
imation techniques of Q to show that 2-local Hamilto- 
nians are sufficient to implement HQC. Furthermore, we 
show that the construction of Q etc., originally devel- 
oped for computational complexity proofs, has a hidden 
gauge freedom and can be viewed as half a holonomic cy- 
cle. This view allows us to improve the adiabatic imple- 
mentation so that we obtain only the desired final state 
instead of the computational history. 

This paper is organized as follows. In section 2, we 
give the direct way to construct an adiabatic equivalent 
of any circuit without encoding the computational his- 
tory. It makes use of only the same number of qubits as 
in the circuit and a running time of the same order as the 
depth of the circuit, as shown in section 3. This in general 
involves n-local Hamiltonian, and section 4 discusses why 
this can be difficult to approximate by studying entangle- 
ment properties of the ground state of Hamiltonians. In 
section 5 we show one way to construct a 2-local Hamil- 
tonian whose worst case run-time scales exponentially 
with n, illustrating a tradeoff between resource require- 
ment and running time. In section 6, we look at how 
2-local constructions can be useful for generalized adia- 
batic algorithms and how computation models that use 
Abelian and non-Abelian geometric phases 0, EJ fall into 
this category. Finally, we return in section 7 to the com- 
putational history approach of Q and note its interesting 
connection to holonomic quantum computing. 

2. A DIRECT MAPPING 

We adopt a general definition of adiabatic computa- 
tion and look for a time-dependent, differentiable Hamil- 
tonian H(s), where < s < 1 is the time parameter, 
such that if (0) is an initial Hamiltonian with a unique, 
easily reachable ground state and if (1) is a Hamiltonian 
with a unique ground state encoding the solution of our 
problem. A quantum circuit can be given in the form 
= U1U1-1...U1 |0), where Ui are unitary operators rep- 
resenting one or two qubit gates. To map this transfor- 
mation into adiabatic evolution, we start with a Hamil- 
tonian ff(0), whose ground state is |0), and we would like 
to have if (s) such that \ip) is the ground state of if (1). 
The most common problem in constructing such an if (s) 
is that the energy gap between the ground state and the 
first excited state varies during the evolution. A small 
gap implies a larger probability for the ground state to 
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be excited, and in turn a longer evolution time if we want 
to compensate for it. 

Our main observation is that it is possible to main- 
tain a constant gap size as long as we keep the Hamil- 
tonian H(s) to be of the form U(s)H(0)W (s). Let us 
be more specific. Suppose the circuit requires us to per- 
form unitary gate U on state |0). Let K = —ilogU and 
U(s) = exp(isK), such that U(0) = 1 and U(l) = U. We 
start with a Hamiltonian H with |0) as its ground state: 

H\n)=E n \n) (1) 
We can add V(s) such that the following is true: 

(H + V(s))U(s) \n) = f(n, s)U(s) \n) (2) 

if 

V{s)U{s) \n) = [U(s),H] \n)+(f(n,s)-E n )U(s) \n) (3) 

This completely specifies V(s), and if /(n, s) = E n , 
V(s) is just, in the original (computational) basis, 
U(s)HU(sY — H. It is clear that as s goes to 1 slowly, 
we obtain U \Q) as our ground state without worrying 
about a shrinking gap. Note that f(n, s) allows us to 
manipulate the gap size. 

Using the idea above, we can now spell out the explicit 
mapping. Given Ui, ...Ui, we first replace the overall time 
parameter s by a series of time step parameters for 
i = Si C [0, 1]. This means: 

#(*) = ( I] U(si))H(0)(l[ W(8i)) (4) 

i=l..l i=l.J 

Let the Hamiltonian at the beginning of the i-th time 

''^ where hj 
^1 



step be H«-V = Y,j hf = J2j hf + h 3 , 

denotes individual local Hamiltonians. J^" and J^" 1- refer 
respectively to terms whose qubits overlap with those 
of Ui and terms that act on different qubits. In this 
notation, we can write V{si) as 

V(s i ) = U i (s i )(^\hf)u}(s i )-^hf ) (5) 

For illustrative purpose, let us consider a typical term, 
where Ui is the controlled-Z gate (which with single-qubit 
gates is universal) acting on the first two qubits, and 
acts on the second qubit as well as some other qubits. 
The matrix representation of hj and Ui for the first two 
qubits looks like 



,(*) 



I hi h 2 

h 3 hi 



\ 



\ 




1 
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\ 


hi h 2 


,Ui = 




h 3 hi ) 






- 1 / 



Then UiisJhfuJisJ-h 



/o 







{e- is *- x )h,2 
\ (e^- 1 )^ J 



(6) 



Note from this example that if H^ 1 ^ = hj is m- 
local, V(si) can be at most (m + l)-local, and this hap- 
pens when exactly one qubit of a two qubit gate Ui over- 
laps with one qubit of h^\ Thus V(s{) can be up to 
rt-local where n is the total number of qubits. We will 
study more closely the complexity and locality of such 
Hamiltonians in section 4. 

Let us look at another specific example using Pauli 
matrices X, Y and Z as basis. Suppose for two qubits we 
start the Hamiltonian: 



H = ZZ - ZI + IZ 



(7) 



where ZI means a Z on the first qubit and identity on 
the second qubit etc. Clearly the ground state is |10). 
With that as our starting point, we can apply a CNOT 
and see how it turns into |11). The recipe above tells us 
the Hamiltonian we need to add is 

V(s) = sin(sir)IY + (1 - cos(sit))IZ 

-sin(sn)YZ - (1 - cos(sir))ZZ (8) 

We can see that as s goes to 1, the new Hamiltonian will 
become H + V(1)=—ZZ + ZI + IZ, whose ground state 
is indeed |11). The IY and YZ terms are zero at the end 
points, as the extra terms in [5| are. 



3. ERROR BOUNDS 

We now check the evolution time required for each 
step according to the Adiabatic Theorem Q- The re- 
sult here is useful for the construction in section 5 and 6 
as well. Under the adiabatic approximation, the evolv- 
ing state is proportional to the instantaneous eigen- 
state of the time-dependent Hamiltonian. Substituting 
this into Schrodinger's equation, this means the time 
derivative does not take one eigenstate to another, i.e. 
(m, Si\ \0, Si) ~ 0, to ^ 0. Thus the correction to 
the approximation must be proportional to (we define 
Ki = —ilogUi below) 

lv-/ , d(U( Sl )H^UHs t )) 
— {m, Si\ — ^r-, |0, s, 



(E m - E )dsi 



1 &m — £jQ 



— {m,Si = 0| 



E m — Eq 
i -<m, «i = 01^10,^ = 0) 



|0, Si = 0) 



(9) 
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where |m, s,) denotes the instantaneous eigenstate with 
eigenvalue E m and Si — U/T. i?^" 1 ) preserves the spec- 
trum of \m, Si = 0), so the contribution to the above term 
is due to Ki . Taking Ui to be controlled-Z as an example 
again, the eigenvalues of Ki are and ir. a(si) is there- 
fore bounded by ir/T. The total time required for the 
step is proportional to the transition probability to other 

states, which according to Q, is bounded by 

for the smallest E m . Remarkably, the error is not only 
independent of total number of qubits n, it is also in- 
dependent of Si, which means further local variation in 
evolution speed is not required to achieve optimal tim- 
ing. Of course, U(si) — exp(isK) is just one arbitrary 
choice we make; there may be other forms of U(si) that 
yield better performance or are easier to implement. We 
should note that it is possible to eliminate the error al- 
together by adding auxiliary terms to the Hamiltonian, 
but this would only be useful for state preparation as it 
generally requires complete knowledge of what we want 
to generate. 



4. LOCALITY OF THE HAMILTONIAN 

In hindsight it should not be surprising that this direct 
mapping yields an n-local Hamiltonian. After all, while 
it is easy to decompose an n-local unitary operator into 
a product of 2-local ones, since 2-qubit gates are uni- 
versal, it is far more difficult to approximate an n-local 
operator with a sum of 2-local operators, even with the 
addition of ancilla qubits. This section is devoted to the 
understanding of this difficulty. 

First we review some results by Haselgrove et al [6j . In 
the authors show how the entanglement of the eigen- 
states of a Hamiltonian is related to its coupling topology. 
Intuitively speaking, if an eigenstate shows strong corre- 
lation between bodies which the Hamiltonian does not 
directly couple (i.e. act nontrivially on all of them as a 
tensor product), the Hamiltonian cannot distinguish very 
well between such a state and other similarly entangled 
states that are orthogonal to it. This results in a small 
energy gap. The following theorem from 6] makes this 
idea concrete and suffices for our purpose. 

Theorem 1 Consider a state and a Hamiltonian 
H whose eigenvalues and eigenstates are Ej and \Ej) 
respectively with j=0..d-l; d is the dimensional of the 
Hilbert space and \Eq) is the ground state of H. Let 
F be the overlap of \ip) with \Eq) and E tot be the dif- 
ference between the maximum and minimum eigenval- 
ues. Then for all density matrices p with eigenvalues 
Pi < P2 < P3--v such that tr(pH) ={ip\ H \tp) , the follow- 
ing inequality holds: 
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The proof is elementary and we will refer the readers 
to the lucid explanation in |(J. Now we may apply this 
theorem on the construction in section 2. 

Claim 2 Let Hq be a 1-local Hamiltonian with unique 
ground state |0). There does not exist general k-local ap- 
proximation for the n-local Hamiltonian UHqW , where 
k < n — 1 and U is a polynomial- sized circuit, such that 
the approximation produces exactly the same ground state 
and first excited state. Specifically, one cannot always 
construct a k-local Hamiltonian which has U |0) as a non- 
degenerate eigenstate. 

Proof: We will start with the case without ancilla 
qubits. Consider the state |i/i) = ^=(|000...) + |lll...)),thc 
n-qubit GHZ state. Consider also a k-local Hamiltonian 
H whose ground state |jEo) = |^>), so F=l. If we choose 

p = |(|000...> (000...| + (11) 

where again indicates n zeroes or ones, it is easy to 
see that tr n —kp=tr n —k \ip) (tp\ for k <n—l, where tr n ^k 
means tracing over any n-k qubits. It then follows that 
for k < n — 1, tr(pH)=(ip\ H \tp) for a k-local Hamilto- 
nian H. Putting this into the inequality in Theorem 1, 
we obtain E\ — Eq=0, meaning that the ground state 
corresponding to Eq is degenerate. 

Now suppose a k-local exact approximation ex- 
ists. Choose the n-qubit polynomial sized circuit, 
V=Hadamard 1 ]JZi~ CNOT hl+1 , i.e. a Hadamard 
gate acts on the first qubit, followed by a series of CNOTs 
on the first and second, the second and the third, and 
so on. Clearly, this circuit acting on the initial state 
|000...) produces Suppose we start with a simple 1- 
local Hamiltonian Hq (e.g. Set Hq = Y^,7=i a z) which has 
|000...) as a non-degenerate ground state. If there exists 
a procedure that exactly approximate UHqII^ with a k- 
local Hamiltonian, k < n — 1, this implies there exists a 
k-local Hamiltonian which has \^p) as a non-degenerate 
ground state. Hence we have arrived at a contradiction. 

To generalize this to the case with ancilla qubits, a 
slight extension of Theorem 1 is needed. Let the ground 



states of the k-local approximation be 



"3h 



where 



"3h 



j = l..m enumerates the degeneracies due to the 
ancilla qubits. This product form is necessary if we 
want the computational qubits to remain as \ip). Let 
P = £j=i P ® \ a j) ( a j \ l m - Setting Eq = 0, it is easy 
to check that tr(p H) = Y^jLi tr{p <£> \aj) (aj \ /m) = 0, 
which would force the ground state degeneracy to be 
2m. This in turn implies that there must be degenerate 
ground states due to states orthogonal to contradict- 
ing the assumption that there is a k-local Hamiltonian 



with 



as the only ground states. 



E )p 3+ i < (1 - F 2 )E t , 



(10) 



So far we have seen that a local Hamiltonian can- 
not have certain states as its ground state, as shown 
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by |g, llfll • This is rather expected as it is well-known 
\U\ that there are quantum states not determined by 
any reduced density matrices. The more interesting con- 
nection we would like to point out here, however, is the 
tradeoff between proximity to a non-local state and the 
energy gap, as apparent in Theorem 1. Since the energy 
gap condition is essential to adiabatic algorithms (while 
some forms of adiabatic theorem without gap condition 
exists, they cannot guarantee the final state to arbitrary 
accuracy |l2j), this places another direct tradeoff between 
accuracy and running-time. Now that an exact approxi- 
mation is not possible, we will look at how close we can 
get. 

In |4| , a 2-local approximation for 3- local Hamiltonians 
is contructed (see section 5). Normalizing the total en- 
ergy to unity, the ground state energy gap for the 2-local 
Hamiltonian scales as 6 3 for a ground state O(S) close to 
the original ground state. In fact we can use Theorem 1 
to make this more precise: If the energy gap scales as 5 3 , 
the ground state for the 2-local Hamiltonian has to be 
at least 0(5 3 ) away from an original GHZ- type ground 
state. This proves that there does not exist an approxi- 
mation scheme better than Q in such a way that the en- 
ergy gap scales, say, logarithmically (i.e. Ofi/logS^ 1 )) 
instead of polynomially with the accuracy 0(5) of the 
ground state. 

Following this idea, we can place some bounds on how 
good the approximation for an n-local Hamiltonian can 
be. For simplicity we will consider a 2-local approxima- 
tion, should one exist, that has a unique ground state, 
ancilla qubits included. Consider the state \(f>): 

|0) = i(|000> + |111>) ® i(|000> + |111)) ® ... (12) 

which is a tensor product of mostly 3-qubit GHZ states. 
It is not difficult to see that there are ~ 2"/ 3 orthogonal 
states to it there are not distinguishable by 2-local terms. 
Thus we can form a density operator of rank ss 2™/ 3 and 
substitute it into the inequality HUfl . This tells us that 
the average energy of these ~ 2™/ 3 states has a gap with 
the ground state that is at most (1 — F 2 ) ~ 0(5). We 
can tighten this bound a little by considering the dis- 
tribution of states. If we start with the 1-local Hamil- 
tonian Hq — X)"=i °z (t ne minimal form required for a 
unique ground state), U H U 1 ' has n eigenvalues Ej = j/n 
with degeneracy n\/[j\(n — Simple counting shows 
that for the lowest ~ 2™/ 3 states, the average energy is 
at least rj nE\/Q. Thus we can tighten the bound to 
E x -E < 0(5 /n). 

Before we conclude this section, we should briefly note 
another line of attack due to IQj. Generalizing be- 
yond the GHZ-type states, states corresponding to non- 
degenerate quantum error correction codes (QECC) also 
turn out to be interesting for the study of local Hamilto- 
nians. These are states with the property that, for some 
constant t usually much smaller than n, any Pauli ma- 



trx operators acting non-trivially on up to t-qubits will 
take the state to a set of orthogonal states. Therefore, 
for a QECC state \x) 1 if the operator H — EI is t-local, 
(H — EI) \x) will be a sum of orthogonal states, imply- 
ing that \\(H — EI) \x) || cannot be close to zero. This 
prevents any QECC states from being even close to any 
eigenstate of a t-local Hamiltonian. Recasting this result 
in our language, we see that there does not exist a k-local 
approximation to arbitrary n-local Hamiltonians UHqU^ 
for sufficiently large n without ancilla qubits , because 
QECC states can also be generated efficiently by quan- 
tum circuits (see references in With ancilla qubits, 
however, cancellation can occur for (H — EI) \x) \a) if \a) 
is not a QECC state. We obtain instead a set of con- 
straint equations that the approximation Hamiltonian 
has to satisfy in order to produce QECC states as an 
eigenstate. 

We hope that the discussion above would be useful for 
further research on not only the possibility of local ap- 
proximation, but also the connection between local prop- 
erties of Hamiltonians and polynomial-sized quantum cir- 
cuits. 



5. A LOCAL APPROXIMATION USING THE 
THREE-QUBIT GADGET 

After an abstract discussion of possible local approxi- 
mations, we will now look at explicitly how an approx- 
imation scheme can be used in an adiabatic algorithm. 
We will use as an example the 3 to 2-local reduction intro- 
duced by referred to as the three-qubit gadget from 
now on. 

Let us begin with some general considerations. In or- 
der to directly map a quantum circuit to an adiabatic 
algorithm gate by gate with some approximation every 
time, we would need L approximations where L is the 
depth of the circuit. To achieve an error within 0(e) for 
the final state, including the error due to adiabatic ap- 
proximation discussed in section 3, we can estimate the 
required accuracy at each step as the following. Consider 
the worst case scenario, when all the errors accumulated 
are in the same direction. We first express the angle be- 
tween the correct final state and the approximate final 
state as 9 = O(yfe) for small 9. The average angle ac- 
cumulated at each step is 9/L because the unitary gates 
preserve angles. The allowed error at each step is there- 
fore 1 — cos 2 (9/L) — 0(e/L 2 ). Hence as long as the 
energy gap size scales polynomially with this allowed er- 
ror, the adiabatic algorithm is efficient. We will see one 
such example in section 6. 

At every step, however, if we repeatedly apply the same 
approximation procedure on the approximate Hamilto- 
nian from the previous step, the energy gap size would 
generally not scale polynomially with the allowed error. 
This is because, as observed in section 4, the energy gap 
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will have to be scaled down by at least a factor of 0(6) 
for an allowed error of 0(6). Repeatedly approximating 
approximate Hamiltonians thus results in an energy gap 
of at most 0(8 L ). This would hold true for any schemes. 

We may now look at the specific scheme based on 4|. 
The authors develop a framework of perturbation the- 
ory that gives sufficient conditions for how one Hamil- 
tonian can approximate another. The basic idea they 
consider is as follows. A 3-local Hamiltonian H 3 can be 
represented as a 2-local Hamiltonian restricted to a cer- 
tain subspace, the intuition being that when the interac- 
tion involves more bodies, we have finer restrictions on 
the eigenspaces. Let this 2-local Hamiltonian be Vi and 
the subspace be S. If we add another 2-local Hamilto- 
nian Hi , such that Hi is zero on S and large everywhere 
else, it is not difficult to see that the lower spectrum of 
Hi = Hi + Vi is close to that of H 3 , as Hi gives penalty 
to states outside of S and restricts Vi to S. 

With this intuition, the next tool we need is a good 
measure of the lower spectrum of Hi. This is provided 
by the self-energy £_ (z) (analogous to the sum of one 
particle irreducible diagrams in field theory) defined as 
follows. First we define the Green function G(z) of Hi 
as 

G(z) = (zI-Hi)- 1 (13) 

Now we define £_(#) by 

G__(z) = (z/_-I]_(z)r 1 (14) 

where G (z) is G(z) restricted to the lower spectrum 

of Hi (not H 2 \). With this definition, proved that 
(Theorem 4, Lemma 9) if 

||£_(*0-.ffe//||<* (15) 

for some operator H e f f , then both the lower eigenvalues 
and the ground states of Hi will be 0(6) close to H e ff. 
From this result, we would have a good approximation 
for H3 if there is an H e f f that is manifestly the same as 
H3 on the computational qubits in the energy range we 
are interested in. 

For any 3-local term H3, Kempe et al propose an Hi 
on ancilla qubits and a Vi coupling the computational 
qubits with ancilla qubits, such that when we calculate 
£_, the above equation is satisfied. This construction 
is called a three-qubit gadget. To apply this to our adi- 
abatic algorithm, we note from section 2 that for each 
2-qubit quantum gate we add to the Hamiltonian, a m- 
local Hamiltonian can become at most m+l-local. This 
means if we start with an 1 or 2-local Hamiltonian and 
apply the three-qubit gadget at every step, we should ar- 
rive at a 2-local Hamiltonian at the end. Let us write out 
the terms explicitly: 

To begin with, the following Hamiltonian on the ancilla 



qubits (playing the role of Hi above) is added: 

fi-3 1 

Hanc = ^X!X! / ®( CT ™l^™2+ CT I Z m l<m3+ Cr l m2 Cr J Z m3- 3/ ) 

i—l m 

(16) 

Terms like of ml are Pauli matrices on ancilla qubits 
identified by three indices: i corresponds to the time step 
which runs from 1 to I; the meaning of the second and 
third indices will become clear shortly. <5 would become 
the error of the 2-local approximation; a smaller 6 would 
correspond to better approximated spectrum and ground 
state. Next we give an stepwise approximation, such that 
given a 2-local Hamiltonian H^ -1 ' at the beginning of 
each time step (see section 2), we find a 2-local pertur- 
bation V'(si) to approximate the possibly 3-local V(si) 
(=UiH^-^Uj - H^-V) when U t is applied. To do this, 
we first write (0 in the following form: 

V( Sl ) = U^s^lhfwlis^-Y^hf 

= Yi — BimiBimiBimS (17) 

m 

where Yi is 2-local and the i?'s are positive semidcfi- 
nite commuting operator acting on three different qubits. 
This decomposition is always possible because the Pauli 
matrix product a a ® a 13 (£> c 7 forms a basis for 3-local 
matrices. If the coefficient of a term is positive, we can 
rewrite the basis term as (1 + <r a ) <g> (1 + <J^) <S> (1 + er 7 ) 
+ 2-local terms; if it is negative, we can use rewrite it as 
-(1 - a a ) ® (1 + a 13 ) (gi (1 + cr 7 ) + 2-local terms. This 
way we arrive at the form of (|17|) . and we can see that 
to is the number of such product terms in the decompo- 
sition. Note that while this decomposition may not be 
obvious in practice, it is a constructive procedure that 
can be done with a classical computer program. Now we 
can construct V'(si): 

V'(si) = K + Y^iS-^B^+Bli + B^) 

m 

- 6~ 2 (B im i <g> af ml + B lm2 ® a x lm i 

+ B im3 ®<jf m3 )} (18) 

where the Pauli matrices in the last sum act on the ancilla 
qubits. Each term in the sum involving three ancillae is a 
three-qubit gadget. In summary, our total Hamiltonian 
is H anc + H(Q) + J^i V( s i)i an( i the error introduced in 
this 2-local approximation at each time step is 0(6). 

To check that £_ satisfies equation (|15fl . put H 2 = 
Hanc, V2 = y'(si), and expand £_ as 

£_(z) = V— + (z - A)- 1 V-+V+- + (z - A)- 2 V-+V ++ V+- 
+ (z - A)- 3 V- + V ++ V ++ V+- + . . . 

where A is the gap of H 2 and V-\ denotes the part of V 2 

that couples the lower spectrum to the upper spectrum 
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etc. We can obtain, after some algebra, 
£_0) = Yi®I anc 

M 

- 6 ^2 Bi m iB im2 Bi m3 (g) (of TOl <g> of m2 ® of m3 ) 

m=l 

+ 0(S). (19) 

Since the i?'s are semi-positive definite, the lowest 
eigenvalue is achieved when (y\ ml ®a^ m2 ®of mZ is replaced 
by 1 (i.e. the ancilla qubits are in (|000) + |111))/a/2. 
H anc has restricted the ancilla qubits to be in the sub- 
space spanned by 1 000) and 1 1 1 1) ) , and we effectively 
recover V(si). We can see that the purpose of those 
tJimi . . . terms is to enforce the product relation among 
the {-Bimi, Bi m 2, Bi m z}. Notice the first excited state is 
also the same as in V(si) because the B terms are positive 
definite, so if the ancilla qubits are in (|000) — |lll))/\/2, 
the increase in energy would be more than that of the 
excited states due to the computational qubits. 

Under what condition will this procedure be ineffi- 
cient? Note that in (|18|l . the original B im i terms are 
multiplied by 6~ 2 . The local reduction requires the ap- 
proximating terms to be very large compared to other 
terms in the Hamiltonian. If reduction is later applied 
repeatedly to terms coupling to ancilla qubits from the 
previous steps, the energy level required for the reduc- 
tion scales exponentially. When we normalize the total 
energy to unity, this equivalently means the gap between 
the ground state and the first excited state shrinks expo- 
nentially. Such repeated approximation could be useful 
when we need to implement a shallow circuit with the 
noise resistant properties of the adiabatic computation 
model and the restriction of 2-local interaction. For ex- 
ample, in conjunction with teleportation circuits, the re- 
peated approximation may not be necessary as we can 
teleport many independently and adiabatically prepared 
unitary operations. For a generally efficient mapping, we 
would need either a procedure that directly reduces an n- 
local Hamiltonian to a 2-local approximation, subjected 
to the constraints described in section 4, or some kind 
of adaptive mapping that exploit structures of specific 
circuits. 

In the next section, we will see how similar repeated 
use of the three-qubit gadget can give rise to an efficient 
adiabatic algorithm once we relax the model. 

6. GENERALIZED ADIABATIC ALGORITHMS 
AND HOLONOMIC QUANTUM COMPUTING 

As mentioned in the introduction, if we are willing to 
relax the property that the Hamiltonian keeps track ex- 
actly what the correct state is, we will have more freedom 
to design the Hamiltonian. Going back to the construc- 
tion in section 1, it is clear that we need different Hamil- 
tonians for the same quantum gate at different stages of 



the computation. Yet we know that the unitary transfor- 
mation due to the application of a time-dependent Hamil- 
tonian over a period of time, U = Texp(i H(t)dt) is 
independent of the state, so why do we need different 
Hamiltonians for different stages? The reason is that we 
have so far ignored the phases of the transformation due 
to the adiabatic evolution. The phases include both dy- 
namical and geometric components: 

(P n (T) = exp(-i [ E n {t)dt + t7n(T)) (20) 
Jo 

where 

J n {T)=i^ (n,t\j t \n,t)dt (21) 

is the geometric phase 0] and E n (t), \n, t) are the eigen- 
values and eigenvectors of H(t). Therefore, if we naively 
apply U(t)H W (t) to a state \tp) without making sure 
that is an eigenstate of H , relative phases can de- 
velop between the eigenstates superposed to form \ip). It 
is not difficult to cancel out the relative dynamical phases 
- all we need to do is to apply — Hq for the same period 
of time (or modulo 2ir). For the geometric phases, they 
are often ignored in an open path evolution since they 
can be gauged away by choosing a different set of basis. 
However, the moment we decompose 1^) into eigenstates 
of H (0), the gauge is fixed; if we choose a different set of 
basis at some other time, we would not obtain U \tp) at 
the end. Therefore the open-path geometric phase must 
be taken into account and cancelled accordingly. This 
gives rise to the following partially adiabatic algorithm: 
To apply two-qubit gate U from a circuit: 

1. Pick a simple, 1-local Hamiltonian Hq, \\Hq\\ < 1. 

2. Apply the \(I + U(t)H U(tf) at any rate from t=0 
to T; [7(0) = I and U(T) = U 

3. Calculate G such that G \n) = j n \n). 

4. Apply \(I- H + 2G/T) for time T. 

This algorithm, of course, does not enjoy property II) 
for all time because of step 4. But if we make good use 
of the transformation due to the geometric phase, such 
that, for example, e lG U(T) = U (which is nontrivial to 
solve since G depends on the path U(t)), the G term 
can be dropped from step 4, and the cancellation can be 
greatly simplified - in fact the cancellation would be the 
same whatever gate we want to implement. 

This is also reminisicent of the Geometric Quantum 
Computation model ||, which uses the Abelian geomet- 
ric phase to implement each gate and requires the can- 
cellation of dynamical phases. The difference is that 
our algorithm uses an open path and thus involves non- 
geometric components. 

Next let us consider how we can avoid having to cancel 
the phases. This is only possible if any state we want to 
apply the quantum gate on is an eigenstate of Hq. But 
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section 4 tells us this state cannot be a unique eigenstate 
without Hq becoming n-local in general, so we will have 
to deal with degenerate states. Without any knowledge 
about the state, we would have to make all 2™ n-qubit 
states degenerate - but this means Hn acts trivially on all 
qubits! The dilemma is solved by adding ancilla qubits - 
we can arrive at a non-trivial U (t)HoU (ty if U(t) couples 
between the computational qubits and the ancilla qubits. 
Notice that C/(0)iJ o C/(0) t = U(T)H U(T) since both are 
trivial on computational qubits, so the Hamiltonian goes 
through a cycle. All the relative phases are accumulated 
between states corresponding to the ancilla qubits and 
do not affect our calculation. We have thus arrived at 
the Holonomic Quantum Computation (HQC) modelQ. 

In order to have a non-trivial control manifold, we can 
use qutrits with states |0) and |1) as the usual qubit states 
and 1 2) for control, or we can add ancilla qubits. In gen- 
eral, identifying the manifold that has the right holonomy 
group and finding the path for each two-qubit transfor- 
mation in the circuit model is very difficult. Recently, 
Tanimura et al 0] settled the mathematical question 
of finding a shortest path given an arbitrary holonomy 
group element in a homogenous bundle, which implies 
that with the addition of just one ancilla qubit, we can 
implement any two-qubit transformation in the space of 
computational states. Let us consider the implementa- 
tion of a CNOT gate as an example of this result. Take 
the Hamiltonian on the ancilla qubit to be 



E 
E x 



We can write the time-dependent Hamiltonian, includ- 
ing the two qubits to be transformed, as 

H(t) = E x e tx ViVle- tx + E e tx V Q Vj e~ tx ; t e [0, 1] 

(22) 



where Vq, V\ are 8x4 matrices: 



I4 denotes the 4x4 identity matrix. We start by preparing 
the ancilla qubit in g round state Eq. To implement a 
CNOT optimally, [lj] found X to be: 



X = 



A B 

-st 



where A and B are: 



.4 



(2 








f 


^ 


2 





B = 


ITT 








1 


1 







-1 


{0 1 








\0 


1 / 



We can see that the Hamiltonian only acts on the two 
computational qubits and the one ancilla qubit, so other 



computational qubits are not affected at all. To do the 
same computation using 2-local Hamiltonians, we can 
now apply the three-qubit gadgets of Q described ear- 
lier. At the end of each cycle, unlike the case in section 
5, the ancilla qubits for the reduction can be discarded 
and reused in the next step. The total number of an- 
cilla qubits required is three times the number of terms 
m Em E>i m iB lm2 B im3 of iJTZl, which is a constant. Fol- 
lowing the same analysis, if we want the final state to 
be accurate up to O(e), the allowed error at each step 
should be 0(e/L 2 ) where L is the depth of the circuit. 
The energy gap required is thus 0(e 3 /L 6 ) and the run- 
ning time is 0(L 12 /e 6 ). Note that this bound may be 
far from tight, and it is quite possible that most circuits 
can be implemented in far shorter time. In any case, we 
have arrived at a fully adiabatic evolution that computes 
efficiently any problem solvable by quantum circuit us- 
ing only 2-local Hamiltonians and a constant number of 
ancilla qubits. 



7. CONNECTION TO THE HISTORY 
APPROACH 

Comparing the above running time estimate to the re- 
sults in 0, 0| , which we will refer to this as the history 
approach, it appears that the latter still seems superior 
even though it produces the entire history of the compu- 
tation instead of only the final output (This just means 
one has to repeat the process O(L) times to get the final 
output or put in identity gates.) It is natural to ask what 
is special about this approach and how it is related to the 
models we have studied. For a circuit of n qubits and L 
gates, the history approach has the following final state: 



L 



inalj — 



UiU, 



l-l- 



|0}<8 \l l L - 1 ) (23) 



where |l'0 i_ ') c denotes the state of L ancilla qubits serv- 
ing as a clock. We refer the readers to 0] for the general 
form of the 3- local Hamiltonian which has the above state 
as its unique ground state. Let us just look at the sim- 
plest example - a circuit with one two-qubit gate U and 
n computational qubits. The initial and final states are: 



\lpin 

^(|0)®|0) C 



= |0>® |0) c ; 
U\0)®\1) C ) 



(24) 
(25) 



where one ancilla qubit suffices for the clock. The corre- 
sponding Hamiltonians are: 



H out = -(/®|0) (0| 
-17® ID (0| c - 



c + i® 

-f7 f 



|1>(1| C 
|1>(1| C 
0) <i| c ) 



(26) 



(27) 
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Actually we have deliberately omitted a term |1) 
|0) (0| c that would make \ipini) the unique ground state 
of H ini . This does not change our analysis as long as we 
prepare the state in |^, n ,) at the beginning. In [3| this 
is the Hin put term, which does not affect the evolution 
equation because it remains zero for all time. Aharonov 
et al shows - for this case, trivially - that slowly interpo- 
lating from H ini to Hf ina i takes \ipini) to \ip final)- Notice 
that the Hamiltonians are highly degenerate in the sense 
that if we replace \ipini) and \^ final) by 

4ni) = V\0)®\0) c ; (28) 
hfinai) = |0> ® |0) c + UV |0) ® |l) c ) (29) 

for some U(n) unitary operator V, the same analysis 
would go through. This gauge freedom is already remi- 
niscent of HQC; let us make the connection more explicit. 
Define the time-dependent Hamiltonian 

ff(t) = i(P(t)/®|0) (0| C + Q(i)/®|1) (l| c 

-R(t)U <g> |1) (0| c - R(t)tf ® |0) (l| c ) (30) 

such that P=R=0, Q=2 gives H ini and P=Q=R=1 gives 
Hfmai ■ It is easy to see that 

\tjj(t)) = {normalization}[a(t) \0)®\0} c +f3(t)U |0)®|l) c ] 

(31) 

is a ground state as long as PQ = R 2 and Pa = R(3. 
Suppose we have completed the original evolution and 
obtained \ipfi na i) (ot, /3 = 1), but decided to keep going. 
We can slowly change from P=Q=R=1 to P=2,Q=R=0 
while maintaining PQ — R 2 (this helps keeping a large 
gap), and we would arrive at U |0) (g) |l) c (a — 0, [3 = 1). 
This is essentially a holonomic cycle if we adiabatically 
rotate the clock qubit back to |0) c . In practice we can just 
relabel the clock qubit if we want to do this repeatedly. 

This gives rise to the interpretation that the history 
approach is in fact half of a holonomic cycle. The cycle 
we have just shown is very similar to the CNOT example 
in section 6, as it requires a 3-local Hamiltonian, which 
too can be reduced by the three-qubit gadget, and one 
ancilla qubit; it is merely a different path in the con- 
trol manifold. Repeated application of the gadget allows 
efficient universal quantum computing, as described in 
section 6. 

We can take this further and interpret the history ap- 
proach for L quantum gates as part of a holonomic cycle 
with a control manifold augmented with the space of L 
ancilla qubits. For the 3-local Hamiltonian case, this is 
not difficult to see. Suppose we have starting Hamil- 
tonian Hi and final Hamiltonian Hf, such that linear 
adiabatic evolution takes starting state 

i>t = |0) ® |0) c (32) 



to 

^ / = — L=, Y, Ui...U Q \Q)®\i) c (33) 

^ L + \=0..L 

where Uq — 1 and we have simplified the notation for the 
clock. These Hamiltonians can be constructed according 
to 0], except we omit the term H input = ^ |1) (1^ ® 
|0) (0| c so that we preserve the gauge freedom mentioned 
above. This does not affect the adiabatic evolution be- 
cause the Hamiltonian does not couple between the de- 
generate spaces. Now we can similarly construct another 
pair of Hamiltonian H' l and H'j corresponding to a re- 
verse circuit U\, W L _ X ... such that it takes 

^ = f7i...£/ |0) ® |0) c (34) 

to 

^ = — L= ]T U L ^...U Q \0)®\i) c (35) 

Simple inspection shows that ij)f and ip'f are the same 
up to relabelling of clock qubits. Therefore, if we imple- 
ment the evolution Hi — > Hf(H'f) — * H[ up to relabelled 
clocks, we can take 

i>i = |0) ® |0) c (36) 

to 

4^ = U L ...U \0)®\0) C (37) 

and thus complete a holonomic implementation of the L 
gates. This opens up the possibility of further optimizing 
the path; a more complete analysis is left for the future. 

DISCUSSION 

We have looked at various forms of adiabatic quantum 
computatation and studied their resource requirements 
as well as possible approximations. One issue we have 
not addressed at all is what noise resistant properties 
different models can have. It would appear that the di- 
rect, non-holonomic approach in section 2 places stronger 
condition on the states, as any deviation requires higher 
energy, whereas the degenerate states in the holonomic 
approach have no protection against transition within 
the degenerate level. This generalization may however 
be simplistic; a better analysis should be with respect 
to particular experimental implementations. We simply 
hope that this paper has provided a more unifying pic- 
ture of adiabatic algorithms that will eventually lead to a 
toolbox experimentalists can refer to for different specific 
applications. 

There are at least a few directions for further studies. 
1) In section 2 we constructed an adiabatic equivalent 
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of any arbitrary circuit. What remain unclear are the 
general properties of such Hamiltonians, namely those 
of the form UHqW where H is a simple Hamiltonian 
and U is a polynomial-sized circuit. If it is n-local, in 
what sense is it simpler than the most general Hamil- 
tonian? Are there things about a circuit that we can 
learn through this corresponding Hamiltonian? 2) We 
discussed in section 4 various constraints on approximat- 
ing UHoW with k-local terms. The general picture is 
still unclear; it would be useful to understand precisely 
under what condition an n-local Hamiltonian can be ap- 
proximated by a 2-local one. 3) We have proposed one 
partially adiabatic algorithm in section 6 and we have 
suggested that the adiabatic construction in section 5 
could be useful as a small section of a larger algorithm. 
It is interesting to investigate what merits, if any, these 
partially adiabatic algorithms possess. 4) For HQC, there 
should be room for improvement. If we implement one 
two-qubit gate at a time, it is important to see how the 
use of three-qubit gadget can be optimized and how a 
tighter bound on the evolution time can be obtained. At 
the same time, the history approach discussed in section 
7 may lead to new class of HQC methods that implement 
many gates together efficiently. 5) In order to build an 
efficient computation model using only 2-local Hamilto- 
nian and adiabatic evolution, we have been naturally led 
to the use of non-Abelian geometric phase. The use of 
holonomy, however, may not be the only option. For 
example, open path non-Abelian geometric phase (ljj 
can be non-trivial as well and may lead to novel ways 
of implementing multiple qubit gates. These questions 
are beyond the scope of this paper, and we believe that 
adiabatic quantum computing remains an exciting area 
to explore. 
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